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phase-space reconstruction 

X. Hachair 1 , R. Braive 1 - 2 , G.L. Lippi 3 ' 4 , D. Elvira 1 , L. Le Gratiet 1 , A. 

Lemaitre 1 , I. Abram 1 , I. Sagnes 1 , I. Robert-Philip 1 and A. Beveratos 1 '* 1 

1 1 CNRS - Laboratoire de Photonique et de Nano structures, 

UPR20, Route de Nozay, 91460 Marcoussis, France 

2 Universite Pans Diderot-Paris 7 

3 Universit de Nice-Sophia Antipolis, Institut Non Lmaire de Nice 

4 CNRS, UMR 6618 

* Corresponding author: alexios.beveratos@lpn.cnrs.fr 

Nanoscale lasers sustain few optical modes so that the fraction of spontaneous emission /3 funnelled 
into the useful (lasing) mode is high (of the order of few 10 _1 ) and the threshold, which traditionally 
corresponds to an abrupt kink in the light in- light out curve, becomes ill-defined. We propose an 
alternative definition of the threshold, based on the dynamical response of the laser, which is valid 
even for = 1 lasers. The laser dynamics is analyzed through a reconstruction of its phase-space 
trajectory for pulsed excitation. Crossing the threshold brings about a change in the shape of the 
trajectory and in the area contained in it. An unambiguous definition of the threshold in terms of 
this change is shown theoretically and illustrated experimentally in a photonic crystal laser. 

PACS numbers: 42.55.Tv, 42.55.Sa, 42.55.Ah 



Recent progress in the fabrication of nanoscale lasers 
has opened the possibility of developing extremely effi- 
cient devices in which even the spontaneously emitted 
photons are preferentially funneled into the lasing mode. 
In these devices, the fraction of spontaneous emission 
going into the laser mode becomes a significant frac- 
tion of unity, in contrast to conventional lasers, in which 
~ 10 -5 . High-/? nanolasers may be obtained in high-Q 
small- volume cavities, such as micropillars [l|, l2| or pho- 
tonic crystal nanocavities [3J-J5|], where the funneling is 
due to the directional enhancement of spontaneous emis- 
sion through the Purcell effect, nanowires [6(, in which 
the high results from the very small number of modes 
they sustain [7[, or surface plasmon lasers (SPASERs) 
[8l4l0l|. where the plasmon field strongly enhances spon- 
taneous emission. High-/3 nanolasers can also be obtained 
in the strong coupling regime, using polariton condensa- 
tion \tj&. 

When approaches unity, the transition from spon- 
taneous to stimulated emission becomes smooth, so that 
the traditional definition of the threshold as an abrupt 
change in the Light-in Light-out (L — L) curve of the de- 
vice is no longer applicable 12]. In fact, for = 1 the 
L — L curve is a perfect straight line, with no break, so 
that from the viewpoint of energy efficiency the = 1 
laser can be considered "thresholdless" [13[ as it exhibits 
the same efficiency for both spontaneous and stimulated 
emission. However, from the viewpoint of dynamical re- 
sponse it is important to be able to distinguish between 
the spontaneous and stimulated regimes, as the large- 
signal direct modulation (gain-switching) bandwidth of 
the laser is different in these two regimes. In addition, 
in high- lasers, thanks to the deterministic channelling 
of spontaneous photons into the lasing mode, stimulated 



emission is initiated very rapidly giving rise to a broad 
direct modulation bandwidth even for large amplitude 
modulation (On-Off Keying - OOK) , a feature desirable 
for high-speed optical data transmission. 

In order to differentiate the spontaneous from the stim- 
ulated regime, it was proposed to monitor the threshold 
through the evolution of coherence, in particular through 
the increase in the first order coherence, as dictated by 
the Schawlow-Townes relation [3j, ll4j (which, however, 
can only be measured for stabilized continuous-wave op- 
eration) or through the decrease of the second order au- 
tocorrelation function, from g( 2 )(0) = 2 below threshold 
(corresponding to chaotic light) to K^lfO) = 1 (corre- 



sponding to coherent emission) [l|, [2j, U, la, llli- How- 
ever both techniques have fundamental limitations. In 
nanoscale lasers, thermal or refractive index fluctuations 
lead to a deviation from the Schawlow-Townes expecta- 
tions, while the transition from g( 2 ) (0) = 2 to unity is not 
abrupt 1] and thus does not permit a more precise deter- 
mination of the threshold than the L — L curve. At the 
same time, the threshold obtainable from the evolution of 
coherence is not related in a simple way to the stimulated 
emission threshold and does not provide direct informa- 
tion on the response dynamics of the laser. In addition, 
a definition of the threshold based on g( 2 ) and the un- 
derlying amplitude fluctuations would be inapplicable to 
lasers displaying non-classical photon statistics, such as 
amplitude squeezed lasers 



15] 



In this Letter, we propose a phase-space reconstruction 
technique that permits an unequivocal identification of 
stimulated emission and provides an unambiguous defi- 
nition of the threshold in nanolasers that can also be used 
for nanolasers with = 1. We illustrate this approach 
through experiments identifying spontaneous and stim- 



ulated emission in the output of a high-/3 (/3 « 0.25) 
nanolaser. 

Our analysis is based on the traditional laser rate equa- 
tions which, for a simple laser, are 16j: 



dN(t) 



= Pf(t) - 7||JV(t) -p 1{l N(t)n(t) (1) 

(// = -r c n(i) + /3 7 ||iV(i) +p 1{l N(t)n(t) (2) 



dn(i) 



where N(i) is the number of excited dipoles, nit) the 
number of photons in the cavity, f3 the ratio of sponta- 
neous photons in the laser mode over the total number of 
spontaneously emitted photons, 7,7 the dipole lifetime, 
r c the cavity loss rate and Pf(t) describes the pumping, 
with the integral of /(£) over time, equal to 1. It should 
be noted that, because of the short round-trip time of 
light in the nanoscale cavity, the cavity decay is fast, so 
that nanolasers are generally characterized by T c > 7|| 
(Class-B lasers). Following [17], we can define iV c = £?- 
the critical number of dipoles necessary to attain the 
threshold. More complicated rate equations, taking into 
account various effects of the gain material, such as ab- 
sorption of ground-state dipoles or non-radiative losses, 
would not alter the general discussion of our approach. 

Our analysis requires that the system be unexcited be- 
fore and after the experiment, in other words that the 
pumping last only for a finite time, so that the number 
of photons and excited dipoles return to zero after the 
pumping is over. To visualize the evolution of the laser 
under Eqs. JTJ) and ([2]), we represent its response in a 
two-dimensional phase-space defined by {dn(t)/dt, n(t)} 
[lji], in a manner analogous to that for the traditional 
representation of transient laser dynamics through the 
parametric curve {N(t),n(t)} in which time is factored 
out [l9(. The advantage of this novel phase-space is 
that, in contrast to N(t) which cannot be measured di- 
rectly, dn(t)/dt like n(t) are both accessible experimen- 
tally, through the signal intensity measured outside the 
cavity, S(t) and dS(t)/dt. As the system returns to its 
initial (ground) state at the end of the excitation pro- 
cess, its trajectory in phase-space is a closed curve. The 
rate equations ([T]) and ^$ are solved numerically, tak- 
ing fit) as a Gaussian pulse much shorter than 7,7 . 
Figure (JTJ,) represents typical phase-space diagrams be- 
low (P = 0.1iV c ) and above (P = 2N C ) threshold, for a 
nanolaser with /3 = 0.1 and T c = IO711. At low pump- 
ing energies, when the laser operates in the spontaneous 
emission regime, the number of photons n(t) and its 
derivative dn(t)/dt both increase, but attain relatively 
small values overall. The phase-space trajectory is thus 
a relatively small loop. On the other hand, at high pump- 
ing energies, when the system operates under full laser 
dynamics (20|, the number of photons n(t) and its deriva- 
tive dn{t)/dt both increase to high values, so that the 
area enclosed by the trajectory increases dramatically. 
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FIG. 1: Phase-space trajectories of the nanolaser response: 
(a) Theoretical trajectories for a nanolaser of j3 = 0.1 and 
r c = 107|| , under pumping below (P = 0.1iV c ) and above 
(P = 2N C ) threshold, (b) Experimental time-traces of the 
output of a nanolaser under impulse excitation, below (P sp = 
1.7/iJ cm -2 ) and above (P s t = 8.6/iJ cm -2 ) threshold, (c) 
Phase-space trajectories corresponding to the experimental 
time-traces above. 



In addition, for Class-B lasers (F c > 711) the shape of 
the trajectory loop changes as soon as stimulated emis- 
sion becomes the dominant dynamical process. Below 
threshold, after n(t) reaches its maximum, the phase 
space trajectory returns to the origin through the neg- 



ative half-plane (dn(t)/dt < 0) quasi-vertically so that 
the overall curve lies mostly in the right-hand half-plane 
defined by dn(t)/dt > 0. In that case — T c n(t) is slightly 
bigger than /3juN(t)(l + n(t)). The number of excited 
dipoles is still large and decays with a characteristic life- 
time 7|i . On the other hand, above threshold, the number 
of excited dipoles drops considerably due to stimulated 
emission and the negative term dominates. Thus the tra- 
jectory makes a long excursion in the negative half-plane 
so that the phase-space curve becomes almost symmetric. 
Hence, even for = 1, the effect of stimulated emission is 
observable in the phase space trajectory and in the area 
it encloses. 

The area enclosed within the phase-space trajectory 
can be calculated through the Gauss-Green formula for 
a parametric curve {x(t),y(t)}: 



A=- 



x(t) 



dy(t) 

dt 



y(t)—\dt, 



which in our case reads, 



.4 



m -*>^ 



dt. 



(4) 



Equation ([4| can be integrated by parts, and since at 
t = and t = oo both n(t) and 



dn(t) 
dt 

pulsed excitation, the area is given by 



.4 



■. m'* 



are zero under 



(5) 



(3) which, using Equation ([2]), can be expressed as 



A = 



(l - N(tj) n 2 (t) + N 2 (t) - 2 (l - N(tj) N(t)n(t) 



dt 



(6) 



where N(t) = N(t)/N c . Below threshold, N(t) < 1 at all 
times, so that the third term in Eq. ([6]) is negative and 
is thus subtracted from the overall phase-space area. On 
the other hand, this term turns positive above threshold 
and its contribution adds to the total area during the 
time it remains positive. Thus, the area enclosed by the 
parametric curve is expected to increase dramatically at 
the transition between spontaneous and stimulated emis- 
sion, even in the case = 1 for which no break in the 
L — L curve can be observed. Figure [2] shows a log-log 
plot of the calculated L — L curve as well as the square 
root of the calculated area enclosed by the phase-space 
trajectory as a function of the pump energy for = {1, 
0.1, 0.01, 0.001} and T c = lO'yn . The square root is used 
to facilitate comparison between the two types of curves, 
since both the signal and the square root of the area 
are proportional to P below threshold. The two types of 
curves coincide at low pumping and display the onset of a 
nonlinear increase (a "knee") at approximately the same 
pumping energies. However, while the L — L curves dis- 
play a level change only for /3 < 1, the phase-space area 
curves display such a change even for = 1, permitting 
identification of stimulated emission in all cases. 

In steady-state lasers the threshold is conventionally 
defined as the x-intercept of the high pumping power 
asymptote in the linear L — L curve, which occurs at 
P = (1 — f3)N c . Clearly, this definition cannot be applied 
to the /3 = 1 laser, as it sets the threshold at P = 0. 
Alternatively, the threshold is often defined as the inflec- 
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FIG. 2: Comparison of the theoretical Light-In Light- 
Out curves (continuous black lines) and the phase-space 
area curves (dashed red lines) under pulsed pumping for a 
nanolaser with F c = rOyy and = {1, 0.1, 0.01, 0.001} 



tion point of the log-log L — L curve, where the number 
of photons in the cavity per dipole lifetime is l/y/0, so 
that stimulated emission dominates the response of the 
laser beyond that point. On the other hand, the onset of 
stimulated emission is associated with the "knee" of the 
log-log L — L curve, where the number of photons starts 



exceeding the value of 1. Thus the '"knee"' of the area 
curve can be defined as the thershold. This definition 
is also valid for (3 — 1 laser and identifies the point be- 
yond which the laser dynamics (including its modulation 
bandwidth) change. 

To illustrate the above considerations, we have inves- 
tigated experimentally the phase-space transition for a 
photonic crystal nanolaser. The cavity consists of a pho- 
tonic crystal double heterostructure 12 ll etched in a 180 
nm-thick suspended GaAs membrane [22], incorporating 
a single layer of self-assembled InAs quantum dots at its 
vertical center plane. Details of the structure and experi- 
mental setup can be found in [23] . The nanolaser is opti- 
cally pumped by a Ti:Sa laser delivering 3 ps-long pulses 
at a 81.8 MHz repetition rate and tuned to 840 nm near 
the energy gap of the wetting layer, in order to reduce 
heating of the membrane. Time traces are obtained by 
a Syncroscan streak camera (Hamamatsu), with a tem- 
poral resolution of 3 ps. Two representative experimen- 
tal time traces are represented on Figure [lib), one m 
the spontaneous emission regime (P sp — 1.7 /iJcm~ 2 ) 
and one in the stimulated regime (P st = 8.6 ^Jcm~ 2 ). 
The decay rates of these two traces permit us to deduce 
7.7 1 = 50 ps and T^ 1 = 13 ps. The phase-space trajec- 
tories were obtained by measuring the cavity signal as 
function of time S(t) and then computing dS(t)/dt af- 
ter applying a numerical low pass filter to smooth out 
the noise in the data. Figure (He) represents the phase- 
space trajectories that correspond to the two experimen- 
tal time-traces of Figure Hfb). The two key features ex- 
pected upon passage through the threshold, namely the 
change in symmetry and the dramatic increase in the 
area are clearly seen. Figure [3] presents a log-log plot of 



the square root of the integrated phase-space area, ap- 
propriately normalized, and compares it with the L — L 
curve. A spontaneous emission fraction of /? « 0.25 can 
be inferred from the intensity ratio between the sponta- 
neous and stimulated emission levels. The "knee" of the 
area plot coincides with that of the L — L curve, giving 
for the threshold Pth = 2.1 /iJcm~ 2 . 

In summary, the emission properties of nanoscale lasers 
with large spontaneous emission coupling factors have 
been studied from a dynamical point of view. In such 
lasers, the transition from spontaneous to stimulated 
emission extends over a wide pump power range in the 
Light in - Light out curve, and the threshold cannot 
be defined easily, especially when j3 — 1. We provide 
a simple and unequivocal identification of laser thresh- 
old through a new indicator: the area A enclosed by the 
phase-space trajectory that corresponds to the impulse 
response of the device. With this definition a threshold 
is defined even when none is visible in the L — L curve, 
as is the case for j3 ~ 1 lasers. This definition is directly 
related to the onset of stimulated emission and thus di- 
rectly pertains to the dynamical response of the laser, 
while at the same time it provides a practical tool for 
quick and easy measurements. 
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FIG. 3: Log-log plot of the Light in - Light out curve (blue 
open circles) and the square root of the area enclosed by the 
phase space trajectory (black circles) as a function of pump 
pulse energy. The dashed lines of slope unity in log-log scale 
are a guide to the eye 
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